Multi-wavelength Observations of PKS 2142-758 during an Active Gamma-Ray State 
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PKS 2142—758 is a flat spectrum radio quasar which emits few, weak but significant 7-ray flares 
in the MeV through GeV energy range. The first flare occured on April 4th, 2010, when the source 
reached a daily flux of (1.1 ± 0.3) x 10" 6 ph cm" 2 s _1 (ATEL #2539) in the 100 MeV to 300 GeV 
range. This flux represented more than an order of magnitude increase over its quiescent flux. 
Since the initial flare, this source has been detected in an elevated state within the same energy 
range from October to November of 2010 and another period ranging from July to August of 2011. 
During the latest flaring period in 2011 a multi wavelength observing campaign was carried out 
using the Ceduna radio telescope, the Australian Telescope Compact Array (ATCA), the TANAMI 
VLBI Array, Swift, the Rapid Eye Mount Telescope (REM), and the Large Area Telescope (LAT) on 
board Fermi. These quasi-simultaneous data were used to construct a broadband SED of this object 
in its rare active state. We present these observations and the resulting SED and some preliminary 
analysis of the constraints they place on the high energy emission from this object. 
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I. INTRODUCTION 

Active galactic nuclei (AGN) are the brightest persis- 
tent sources of electromagnetic radiation, from low en- 
ergy radio waves to high energy 7-rays. It is believed 
that AGN are powered by accretion of matter onto a su- 
permassive black hole [9 j. AGN can produce large scale 
jets of plasma that can extend for thousands of parsecs. 
Blazars are the most luminous and violently variable sub- 
class of AGN whose relativistic jets are believed to be 
pointed close to our line of sight [T5] . Blazars are also 
the most common type of AGN to be detected at GeV 
energies. The method by which blazars produce 7-rays 
is uncertain and there are a number of different possible 
explanations. 

Models of high energy blazar emission can be divided 
into two major subgroups, leptonic and hadronic models 
[3]. Leptonic models focus on a source of seed photons 
being Compton upscattered to higher energy levels by 
relativistic electrons in the jet. Leptonic models can be 
further subdivided by the source of seed photons, one 
candidate source is synchrotron photons produced by jet 
electrons themselves (synchrotron self Compton scatter- 
ing). Another possibility is that the seed photons come 
from a source external to the jet (external Compton scat- 
tering). The accretion disk itself, the broad line region 
and the dusty torus are all possible examples of photon 
sources for external Compton scattering. Seed photon 
sources do not necessarily have to be confined to the host 
galaxy, cosmic microwave background photons can also 
be up scattered to higher energies [4]. Hadronic models 
assume that a significant fraction of the jet power is con- 
verted into acceleration of protons. These protons reach 
the threshold for proton pion production and synchrotron 
supported pair cascades will develop [TO] . 

Quasi simultaneous data on blazar emission at multiple 
wavelengths across the electromagnetic spectrum will al- 
low us to distinguish between the different models of high 
energy production within blazars [3 . The launch of the 
Fermi 7-ray space telescope has created new opportuni- 
ties for studying blazars at high energies. The LAT is 
more sensitive and has better energy and spatial resolu- 
tion than any previous telescope in the MeV-GeV energy 
range. Fermi has put better constraints on blazar vari- 
ability within its observing window and has shown that 
it is not unusual for blazars to exhibit fast variability on 
timescales of days or even hours at GeV energies. Fermi 
operates in an all sky scanning mode which makes it an 
ideal rallying point for staging multi instrument observ- 
ing campaigns on sources in response to their changing 
7-ray state. 

The brightest AGN in Fermi 's observing window are 
studied in detail using a broad band multi wavelength 
approach (e.g. [I], [H]). However relatively little atten- 
tion has been given to blazars that are not the most ex- 
treme 7-ray emitters. PKS 2142-758 is a flat spectrum 
radio quasar at a redshift of 1.139 [8 and it is not among 
the brightest Fermi detected extragalactic sources. Typ- 



ical 7-ray flares of this source reach a flux of around 
1.0 x 10~ 6 ph cm -2 s _1 in the 100 MeV-300 GeV energy 
range and they happen infrequently. This object will be 
part of a broader study of blazars which are not the most 
extreme 7-ray sources. The aim of this study is to inves- 
tigate why some blazars can be very 7-ray bright while 
other are not despite their similarities at other wave- 
lengths. 



II. 7-RAY PROPERTIES 

All the 7-ray data presented in this paper were 
recorded by the LAT. The LAT is a high energy 7-ray 
telescope with a wide field of view designed to cover 
the energy band from 20 MeV to greater than 300 GeV. 
The LAT primarily operates in an all sky scanning mode 
sweeping its field of view over the entire sky every three 
hours. Likelihood analysis is used in order to determine 
the gamma properties of any detected source for any time 
range during the LAT's mission lifetime [2 . All LAT 
data analysis was unbinned and done using the 09-24-00 
version of the science tools and the P7SOURCE_V6 IRF. 

The light curve (Figure 1) shows the 7-ray flux of PKS 
2142-758 in the 100 MeV to 300 GeV energy range aver- 
aged into weekly time bins. A circular region of interest of 
10 degrees was used for the analysis. The light curve be- 
gins on August 4, 2008 and ends on September 28, 2011. 
Upper limits on the source's flux are calculated when the 
test statistic (TS) [TT] for that week is below 9 (9 roughly 
corresponds to a 3 sigma detection of the source). The 
three flares originating from this source peak at around 
1 x 10 -6 ph cm^s- 1 in the 100 MeV to 300 GeV range. 
Our multi wavelength campaign was carried out during 
the last flaring period in July and August of 2011. 

The LAT spectrum (Figure 2) was made using a cir- 
cular ROI of 10 degrees and data from July 21, 2011 to 
August 18, 2011 which corresponds to a time period when 
the source's daily TS across the 100 MeV to 300 GeV en- 
ergy band was greater than 25. Initially a likelihood anal- 
ysis was run over the whole spectrum to get a power law 
fit for the 7-ray emission. The data were then divided 
into 13 logarithmic energy bins. Using the initial power 
law as a starting point for the analysis the flux was de- 
termined in each of the energy bins. The horizontal bars 
on each of the data points represent the width of the 
energy bin and the vertical bars are one sigma errors. 
Points with arrows pointing below them are upper limits. 
Upper limits on the target source's flux were calculated 
when the test statistic within the energy bin was below 
25. 



III. SWIFT OBSERVATIONS 

Swift is a multi wavelength space-based observatory 
that has three different instruments on board, the Burst 
Alert Telescope (BAT), the Ultra Violet Optical Tele- 
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FIG. 1: 7-ray light curve of PKS 2142-758 1 week time bins and test statistic < 9 for upper limit threshold 
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FIG. 2: LAT energy spectrum, the model displayed is the 
power law fit for this source across the entire 100 MeV to 
300 GeV energy range 



scope (UVOT) and the X-Ray Telescope (XRT). Swift 
uses momentum wheels to rapidly change the direction 
in which it is pointing [7]. The observatory's primary 
objective is to observe 7-ray bursts making its ability to 
slew rapidly a necessity, however this ability has become 
extremely useful for studying blazars due to their rapid 
variability at high energy. 

Swift observed this source with the UVOT and XRT 
on August 10, 16, 18, and 23 2011 with exposure 
times of 2085 s, 4890 s, 3666 s, and 1441 s respec- 
tively. The most current versions of the Swift data re- 
duction software (07Jun2011_V6.11) and calibration data 
(caldb.indx201 10725) were used. This resulted in four 
UVOT observations with 3 different filters J7, UVW1 and 
UVW2 which correspond to wavelengths of 3465, 2600 
and 1928 A respectively. Aperture photometry was done 
using the standard UVOT data reduction task UVOT- 
SOURCE. A circular source region of radius 2.5" and 
background region of radius 12.5" were used for the cali- 



bration. UVOTSOURCE also performs coincidence loss 
corrections on the data. Magnitudes have been corrected 
for Galactic extinction using the method described in [6] . 

The XRT observations were made in photon count- 
ing mode and a 23.6' x 23.6' image is available for each 
pointing. The source is clearly visible with no other 
sources in the region of interest for each pointing. Identi- 
cal circular source and background extraction regions of 
60" and 119", respectively, were used for all four point- 
ings. Source count rates of 0.054(5) counts/s, 0.040(3) 
counts/s, 0.049(4) counts/s, and 0.043(6) counts/s were 
measured. In order to increase the signal to noise ra- 
tios for spectral modeling the spectra of the individ- 
ual observations were averaged into a single data set. 
An overall response file was created weighing the aux- 
iliary response files by photon number and using the ap- 
propriate response matrix from the calibration archive, 
swxpc0tol2s6_20010101v013.rmf. The averaged spec- 
trum was rebinned by combining every 25 channels below 
3.5 keV and every 50 channels above. The 12 ks spec- 
trum was then modeled in the energy range from 0.3 eV 
to 8keV with an absorbed power law. This resulted in 
a good fit with best fit parameters of Nu — 6.6 x 10 -20 
cm -2 and Gamma=1.50^o.*20 f° r the. neutral hydrogen 
column density and the power law index, respectively 
(error confidence level 90%) and a reduced x 2 value of 
1.2 for 10 degrees of freedom. 



IV. NEAR INFRARED OBSERVATIONS WITH 
REM 

The Rapid Eye Mount (REM) telescope is a 60 cm tele- 
scope observing at infrared and optical wavebands. It is 
able to react quickly to changes in the sky. It is located 
on the La Silla premises of the ESO Chilean Observatory 
making it ideal for observing highly variable southern 
sources. Observations of PKS 2142-758 were made on 
August 10, 2011 shortly after a flare in gamma-ray emis- 
sion, with 150 s exposure time. To calibrate the data, 10 
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reference stars from the 2MASS catalogue were used with 
J magnitudes ranging from 11.10 to 14.9. The linear fit 
between instrumental and nominal magnitudes was al- 
ways good, with a slope of nearly 1, and a dispersion 
of about 0.07 magnitude. Two near IR flux measure- 
ments of PKS 2142—758 were made resulting in magni- 
tudes of 15.05 ± 0.16 in the J band and 14.92 ± 0.31 
in the H band. The errors reported are computed by 
IRAF/apphot and they are based on the Poisson statis- 
tics of the net counts and the detector noise: this method 
generally underestimates the actual error in the source's 
magnitude. The source was rather faint for REM, so that 
the total net counts were rather low and the strongly 
structured background light leads to an additional er- 
ror of 25-30 % (0.3 magnitude) in the flux measurement. 
Our total error for the source flux is about 50% at both 
frequencies. 



V. RADIO OBSERVATIONS 

The radio data reported here were obtained from two 
different monitoring programs one using the Australia 
Telescope Compact Array (ATCA) and the other using 
the Ceduna telescope. 

ATCA has been monitoring TANAMI sources at a 
number of radio frequencies. ATCA data were reduced 
with the Miriad software package, following standard 
practices for continuum data. Flux densities were mea- 
sured using the triple amplitude quantity, and are cali- 
brated against B1934— 638 for frequencies below 25 GHz, 
and against Uranus for higher frequencies. One flux den- 
sity measurement was made for each 2 GHz band, and 
represents the flux at the centre frequency, assuming that 
the spectral index across the band does not vary. For this 
study, PKS 2142—75 was observed on August 30th at 5.5, 
17, 19, 38, 40 GHz. 

The University of Tasmania operates a 30 m radio tele- 
scope at Ceduna in South Australia (e.g. [32]). This tele- 
scope has been monitoring TANAMI sources at 6.7 GHz 
typically once every two weeks. The data are collected 
in four scan blocks, two scans in right ascension and two 
in declination. A block is rejected if any of its scans fail 
a Gaussian fit. The data are corrected for gain-elevation 
effects, and a pointing correction applied. B1934— 638 is 
used as the primary flux calibrator. Variation in system 
temperature due to ambient temperature changes over 
the course of 24 hours necessitates daily averaging of the 
data. The data used in the construction of this SED were 
observed on July 30th, 2011. 



VI. SED MODELING 

The data from each instrument are plotted on the spec- 
tral energy distribution (SED) in Figure 3. An initial at- 
tempt at leptonic SED modeling was performed accord- 
ing to the method described in [5 j. This model assumes 




v[Hz] 



FIG. 3: Quasi simultaneous SED of PKS 2142-758 while in a 
7- ray active state. 

that the radio flux is non-thermal synchrotron radiation 
emitted by electrons which are isotropically oriented in 
the co- moving frame of the jet in a randomly oriented 
magnetic field in the jet plasma. Optical and UV emis- 
sion is modeled as thermal emission from the accretion 
disk. The X-ray through 7-ray emission is produced by 
Compton scattering of seed photons produced in the dust 
torus. The contribution of the synchrotron self Compton 
effect is very small and plays almost no role, so external 
Compton processes are the primary source of high en- 
ergy emission for this source. The variability timescale 
derived from the LAT light curve constrains the size of 
the 7-ray emitting region. The synchrotron self absorp- 
tion process is included in this fit so the small region that 
is consistent with the variability timescale cannot explain 
the radio points. Table 1 lists the parameters that went 
into the model. 



VII. CONCLUSIONS AND FUTURE WORK 

The next step in this project is to construct and model 
an SED of the source when it is in a quiescent state. Com- 
parison of the fitted parameters in the active and qui- 
escent states should provide constraints on the gamma- 
ray production mechanism in this source. Ongoing VLBI 
(Very Long Baseline Interferometry) monitoring by the 
TANAMI program [13] will allow us to resolve the jet 
structure of this source and follow its kinematics. In par- 
ticular, apparent proper motion measurements will pro- 
vide a much stronger constraint on the Doppler factor 
than is possible with the SED modeling. 
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TABLE I: SEP Model Parameters 
Variability Timescale 5 x 10 5 s 
Radius of Emission Region = 2.8 x 10 17 cm 
Doppler Factor = 40 

Magnetic Field = .23 G 

Broken power law electron distribution: 

Minimum electron Lorentz factor = 2.3 
Break in Lorentz factor = 230 
Maximum electron Lorentz factor = 9 x 10 3 
Low energy electron spectral index = 2.0 
High energy electron spectral index = 3.2 
Dust torus: 

Dust temperature = 1 x 10 3 K 

Dust Luminosity = 2 x 10 45 erg/s 

Radius of Dust Torus =16.2 parsecs 

Accretion disk emission: 

Black Hole Mass = 1 x 10 9 solar masses 

Accretion Rate =1/12 eddington rate 

Inner Radius = 6R 9 (R g = 1.5 x 10 14 cm) 

Outer Radius = 1 x 10 4 R^ 

Luminosity of the disk = 6.5 x 10 45 erg/s 

Equipartition: 

Jet power in magnetic field = 5.0 x 10 46 erg/s 
Jet power in particles = 2.9 x 10 45 erg/s 
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